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摘 要 :对 瞬 态 扭转 激 振 荷载 作用 下 层 状 饱和 土 中 端 承 管 桩 的 动力 响应 问题 进行 了 研究 。 首 先 借 
助 Laplace 交换 和 分 离 变 量 法 解 夺 Biot 波动 方程 , 土 层 层 间 动力 相互 作用 简化 为 分 布 式 线 弹 簧 ,并 
结合 柱 - 土 交界 面 上 的 应 力 和 位 移 连续 边界 条 件 , 导 出 了 各 节 段 管 柱 在 Laplace 变换 域内 扭转 动力 
响应 的 解析 表达 式 ; 然 后 ,根据 阻抗 函数 传递 原理 ,并 对 得 到 的 解 进行 Laplace 数值 逆 变 换 , 最 终 获 
得 了 瞬 态 扭转 荷载 作用 下 层 状 饱和 土 中 管 柱 动力 响应 的 时 域 解答 。 研 究 表 明 : 桩 顶 作用 三 角形 或 
< 十 半 正 弦 荷 载 情况 下 ,最 大 桩 顶 捏 转角 产生 时 刻 落后 于 最 大 荷载 对 应 时 刻 ; 且 当 桩 周 土 渗透 系数 处 于 
) 相 对 较 高 和 较 低 范围 内 , 桩 顶 转 角 几 乎 保持 不 变 , 在 中 间 范 围 内 则 随 渗透 系数 增 大 而 增 大 ; 硬 夹 层 
a 的 存在 将 使 柱 身上 部 扭矩 增 大 ,下 部 扭矩 减 小 , 软 夹 层 下 规律 相反 。 
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Abstract : Dynamic response of an end-bearing pipe pile subjected to transient torsional loading in layered 
saturated soil is theoretically studied in this paper. First, Biot's wave equation is decoupled by recourse to 
Laplace transform technique and separation of variables,and a series of distributed linear springs are uti- 
lized to simulate the dynamic interaction. between the soil layers. By combined with the continuity and 
boundary conditions of pile-soil interface, the torsional dynamic response of an arbitrary pile segment is 
further derived in the transformed domain. Then , according to the principle of the recursion of impedance 
functions and by means of inverse Laplace transform , the time-domain solutions of pipe pile under transient 
torsional loading embedded in a layered saturated soil are finally obtained. Numerical results show that the 


time corresponding to maximum twist angle atop the pile lags behind that of maximum loading when sub- 
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jected to triangular or semi-sinusoidal loads. The twist angle atop the pile almost remains constant in the 


case of outer soil with relatively high or low permeability coefficients , whereas it gradually increases in the 


middle range of the permeability. The hard interlayer soil increases the torque of upper part of the pile and 


decreases that of the lower part, while the effect of soft interlayer soil is inverse. 


Key words :pipe pile;layered saturated soil; torsion ; transient response 


在 地 震 、 风 荷载 ,波浪 荷载 和 船舶 撞击 等 冲击 作 
用 下 , 当 管 桩 基础 不 对 称 或 受 偏心 水 平 荷载 时 ,将 产 
生 扭 转变 形 。 此 外 ,即使 对 于 结构 对 称 的 桩 基础 ,在 
诸如 斜 和 射 SH 波 引起 的 地 面 转动 分 量 影响 下 也 将 
发 生 扭 转 振动 。 因 此 ,研究 瞬 态 扭转 荷载 作用 下 管 
桩 动力 响应 特性 ,对 近海 平台 \ 高 层 建筑 和 高 压 输电 
塔 等 桂 承 结构 设计 具有 重要 意义 。 

目前 ,国内 外 学 者 借助 理论 .试验 和 数值 模拟 等 

多 种 方法 对 桩 的 扭转 振动 进行 研究 ,涉及 瞬 态 和 稳 
态 两 种 动 荷 载 形 式 , 单 层 均 质 弹性 和 多 层 非 均 质 铬 
弹性 等 圭 体 特性 ,但 主要 限于 实心 桩 "”) 。 而 对 于 管 
枕 革 础 ,尤其 是 大 直径 管 枉 ,要 同时 考虑 柱 世 十 和 柱 
局 号 的 共同 作用 。 文献 [3] 研 究 了 桩 顶 简 谐 荷载 作 
用 了 B 均 质 线 符 弹性 土 中 管 桩 的 扭转 动力 响应 。 
文献 [4] 引 入 附加 质量 模型 ,探究 士 塞 惯性 .阻尼 效 
应 条 土 塞 与 管 柱 之 间 位 移 相位 差 对 管 柱 扭 转 振动 的 
影响 。 而 对 于 饱和 土 中 管 桩 - 土 动 力 相 互 作用 ， 
文献 [5] 采 用 Biot 饱和 多 孔 介 质 波 动 方 程 ,导出 三 
维 轴 对 称 条 件 下 均匀 各 向 同性 饱和 土 中 端 承 管 桩 受 
谐 和 扭转 荷载 时 的 解析 解 。 文 献 [6] 研究 了 横 观 各 
向 加 性 饱和 土 中 PCC 管 桩 的 受 扭 振动 问题 。 
斧 然 而 ,地 基 士 在 天 然 沉积 作用 和 人 为 因素 影响 
下 6 往往 具有 明显 的 成 层 性 ,导致 土 体 物理 力学 特性 
沿 深度 表现 出 非 均匀 性 。 因 此 ,将 桩 周 土 体 视 为 层 
状 介质 更 符合 工程 实际 。 关 于 成 层 土 中 桩 动力 响应 
问题 ,文献 [7] 基 于 Novak 薄 层 法 原理 ,应 用 传递 矩 
阵 , 研 究 了 层 状 土 中 多 缺陷 柱 在 瞬 态 扭转 荷载 下 的 
动力 响应 。 文 献 [8] 基于 虚拟 的 土 桩 模型 ,对 成 层 
土 中 管 桩 稳 态 扭转 动力 特性 进行 研究 。 文 献 [9] 导 
出 了 成 层 黏 弹性 土 中 管 桩 桩 顶 复 刚度 的 解析 式 。 文 
献 [10] 利用 Winkler 模型 模拟 土 体 沿 深 度 的 分 层 
性 ,并 将 内 外 土 离散 为 多 个 环形 区 域 ,对 番 弹 性 土 中 
管 柱 的 扭转 动力 响应 进行 研究 。 文 献 [11] 考虑 桩 
周 土 三 维 波动 效应 ,推导 了 成 层 亿 和 土 中 桩 的 竖 向 
振动 阻抗 。 基 于 以 上 研究 , 稳 态 荷载 作用 下 成 层 土 
中 桩 土 动力 相互 作用 的 研究 已 较为 深入 ,但 层 状 饱 
和 土 条 件 下 管 桩 瞬 态 扭转 动力 特性 还 尚未 研究 。 


本 研究 基于 Biot 多 孔 介 质 波 动 方程 ,考虑 成 层 
地 基 土 波动 效应 ,以 及 现 浇 混凝土 管 桩 在 成 桩 过 程 
中 可 能 出 现 的 缩 径 和 扩 径 现象 ,通过 分 布 式 线 性 弹 
簧 模拟 各 土 层 间 相 互 作用 ,推导 了 瞬 态 扭转 荷载 作 
用 下 层 状 饱和 土 中 端 承 管 桩 动力 响应 的 解析 表达 
式 , 最 后 利用 Laplace 逆 变 换 求 得 了 桩 身 扭转 角 和 所 
和 矩 的 时 域 响应 ,并 分 析 了 和 荷载 作用 时 间 V Ja fi 
性 和 渗透 系数 等 对 管 桩 动力 响应 的 影响 。 


1 基本 假定 与 方程 建立 


1.1 计算 模型 与 基本 假定 


图 1 为 瞬 态 扭转 荷载 作用 下 层 状 饱和 土 中 管 桩 
的 计算 模型 。 其 中 桩 端 髓 固 , 任 意 激 振 扭 矩 m OD) TE 
用 于 桩 项 , 桩 长 、 土 体 厚度 均 为 有 。 土 体 分 为 N jz, 
各 层 厚 度 分 别 为 请 hs his ,hy; 相 应 地 ,将 管 桩 自 
下 而 上 划分 为 编号 1,2,i,… ,NN 的 入 个 节 段 。r; ,和 
ri 分 别 代 表 第 i (1 iN) 桩 段 的 内 径 和 外 径 ( 下 
标 2 对 应 桩 蕊 土 ,1 对 应 桩 周 土 ,下 同 )。fi,1(z,1)、 
fi.2(z,t) 分 别 代 表 桩 和 映 单位 面积 上 受到 第 i JE ETSI 
土 和 桩 忌 土 的 环 向 切 应 力 。 


刚性 支承 
图 1 有 瞬 态 扭转 丛 载 作用 下 层 状 饱和 土 中 管 桩 计算 模型 
Fig.1 Computational model for a pipe pile embedded in 


layered saturated soil under transient torsional loading 


本 研究 利用 分 布 式 线 弹 簧 模拟 各 土 层 间 相 互 作 
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用 ,将 第 i+1 土屋 和 i-1 土 层 对 第 i 层 土 的 作用 分 
别 简化 为 刚度 为 5 和 ka HAE, E kai 
此 外 ,本 研究 在 建立 方程 时 采取 如 下 假定 。 

1) 假设 土 颗粒 不 可 压缩 , 桩 周 ( 芯 ) 土 为 各 向 同 
性 的 均匀 饱和 多 和 孔 弹 性 介质 。 

2) 忽 略 土 体 径 向 和 竖 向 位 移 , 假 设 只 发 生 环 向 
位 移 。 

3 ) 桩 端 为 刚性 支承 , 桩 周 土 表面 为 自由 边界 ,无 
前 应 力 \ 正 应 力 。 

4) 振 动 为 小 变形 , 管 桩 与 桩 周 ( 芯 ) 土 之 间 完 全 
连续 接触 , 即 桩 土 交界 面 两 侧 应 力 、 位 移 连 续 。 


1.2 方程 建立 
对 任 一 土 层 i(1<i<NN) , 令 土 骨架 ( 固 相 ) 的 环 
[PIRE uu Cr z 0) , 液 相 相 对 于 固 相 的 环 向 位 移 为 
we ,z,1) Hp Biot 波动 方程 , 土 体 扭转 振动 控制 方 
FE 
nalr z De 


GC Ouz) - Gus 
2 


= dein 


si,j 


ð Pins z,i) o wg (r,z,t) 
CA P ag +p: ag (1) 
ÅP MERE p= (1 — n) p, + noz pi HIRKA 
EBP ALDRE, n HIRR; Cu G = 1,2 分 别 
APB SER EAE : 土 ) 为 土 的 前 切 模 量 ; rz 分 别 为 
(BRURIAE IL d cR y? si n no 


ror oz 
第 i 层 土 的 渗流 连续 方程 可 表示 为 
Ee (re) 


9 Prou) pr à wg (r,z,t) 
Lf) +p; 2 Ret 2 
| ot ðt n ðt 


=0 


(2) 
ui^ k, IS Darcy 渗透 系数 ;g 为 重力 加 速度 。 
i 段 桩 扭转 振动 方程 为 
n t) 


Gala 2 T -2mr afia st) 一 -2mr aia (2,0) = 


0 ^e;z,t) 
pul ot (3) 


5C, 6, uio; C2, 0 和 p, 分 别 表 示 第 i ER 
量 、 极 惯性 矩 .扭转 角 和 密度 , 且 几 = mra 702272. 
1.3 边界 条 件 

建立 局 部 坐标 系 来 描述 桩 - 土 边界 条 件 ,原点 设 
在 第 i BEBE (i =1,2,…, N) 桩 顶 , 桩 底 坐 标 设 
为 h,, 则 
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1) 桩 周 土 水 平 无 限 远 处 位 移 应力 趋 近 于 零 , 即 
Tgi i (r—00,2,t) 2 ug í (r—00,2,t) =0 (4) 
2) ies E r 20 处 位 移 为 有 限 值 , 即 
luga (r =0,z,t) |<% (5) 
3) 土 层 顶部 边界 条 件 为 
Qus (r,z,t) -Faiy vE 
7 F "i UP emu 
SN-10154z21 i 
Tay; (0,27 0,1) 0 (j 21,2) (6) 
4) 土 层 底部 边界 条 件 为 
ug ;(r,z zh, ,t) -0 (j=1,2), 
Ug ;(r,z,t) Ev T B 
rH RA 
(122,3, Nu 21,2) (7) 
5) 桩 顶 与 桩 底 边 界 条 件 为 
Aen] mm) uu 
àz Pu sa (i21,2,-,N) (8) 
Qi(z zh, ,t) -0, 
SXODNEACOLS MEM 
| e | 本 仁政 生前 
(9) 
式 中 :5 表示 第 ; 桩 段 底 部 的 反 力 系数 ;m; (1) 为 第 
i+1 桩 段 对 第 i 桩 段 顶 部 作用 的 扭矩 。 
6) 桩 土 界面 位 移 及 应 力 连续 条 件 为 
Ugi (r-rii,2,0) = PiCz,t)T a (10) 
Ugia (T =T; 2,2,6) = pi(2,t)T;, (11) 
fia (2,t) = -Tai (TET ,2,1) 
Oug(r,z,t) ugi(T,z,t) | 
= =G : — + 
a ðr r | zii 
(12) 


fiat) = Tgi2(T 三 六 2 TM) 
=G [an teat) 
CN EL ism 
2 


or r 


3) 
2 问题 求解 


2.1 层 状 饱和 土 振动 方程 求解 


式 (1) 和 式 (2) 两 边 分 别 进行 Laplace 变换 ,得 
a lo "d U,(r,z,s) E 
Cs p Tum gz JU, LE ee ap 


ps PLC rias) t pis ^ Wa(r,z,5) (14) 


PE Wa(r,z,s) +p Us(r,z,s) + Pr (7,z,s) =0 
n 


[m 
(15) 
式 中 :Us(r,z,s) War z,5) ERI ug Cr,z, t) 和 
war ,z,t) 的 Laplace 变换 式 ;s 表示 变换 域内 频率 。 
式 (15) 代 入 式 (14) ,消去 Wu(r,z,s) 可 得 


U.. 2 pen 
V^U,(r,z,s) 加 ct) zs I b ; 
r 


Us(r,z,s) (16) 

式 中 ,2Z,(z) = C;sin(a;z) + D;cos( oz) o 
采用 分 离 变量 法 , 令 Us (7,z,s) 2 ZG) RIT), 

代入 式 (16) ,可 得 

Z;(z) 2 C;sin(ajz) + D;cos( az) (17) 
R;(r) 2 AK, (q;r) +B (qur) (18) 
WRO K(，) 分 别 为 第 一 类 和 第 二 类 一 阶 变 
JcMessel BRA oq; A; B; C; D; HOS CERE CC FH 


boi -— "m pisn 
Q= + 一 一 o 
71e Ga B 
t | Pr 


C9 由 此 可 得 土 的 环 向 位 移 为 

CB CER) 7 [A4;;Ki(q;r) € Bili(qir)] * 

A [Ciysin(æ, jz) +D jcosCaiyz)] — (19) 

CAI 分 别 对 土 层 边界 条 件 式 (4) ~R (7) 进行 
Lablace 变 换 。 当 i=1 时 ,将 式 (19) 代 入 Laplace 变 
换 后 的 式 (4) ` 式 (6) 第 一 式 以 及 式 (7) 第 一 式 ,可 得 


f... G, 
By 0,H tan ( ai, hi) 2 xi EIC -1,2, 
È stl,l 
3,29) ,进一步 可 得 第 一 层 桩 周 土 环 向 位 移 为 
Used rns) = Y, Asia Cmar) sinl æn (2 = h, )1 
m=1 


(20) 

同 理 , 当 i USELESS EI 时 ,可 得 有 =0; 当 
i=N 时 ,可 得 By =0, 则 

Ugal rizs) = > A a miar) sin( einz F Pmi) , 


m=1 
(i e2,3,9 1) (21) 


Ui (7,2,s) 三 > Anr aKı (dny ar) cos( aniz) (22) 
meal 

式 中 

Oi i (kaii + kaii) Gai 

(Ami Gaa) — kaiiEsia 


PPs) p, 
az Gi 


h;) = 


tan (o5, , (m x 1,25 


2 
2mri, 


Gi Ju m= 


s d, (z,s) = 


m= 


Gai > Amim aKa diti )sin[ dus Ct -= h,) ] + 
1 


3i 52,3, N-1), 
kana 
tan(&ny hy) = (m=1,2,3,.…) » 
Anyi Cw. 


Pmi,l = arctan( Oi Cu/ aia ) 
其 中 m: Aa dvwi 为 一 系列 待定 系数 。 
同 理 , 桩 芯 土 环 向 位 移 为 


Ug 2(r,2,s) z > Baal ldm or)sin[l anm a(z- h)] 
mal 
(23) 


oo 
U, 2 ( r,z,8) m £ B, 4l (quiar) sin( Qi 22 T Q2) , 


m=1 


(i =2,3,.…,N-1) (24) 


ie 
Uv. ( 7 和 ,5) 三 > Bal, ( Qual) cos( Ow 22) 
1 


(25) 
式 中 
Qm 26» 
ku 
Q mi2 ( koia + kh: 2) Gs 
( Q mi,2 Gaa) a kai okaio 
(m=1,2,3,.…), 


tan( à, 5h) = = (m=1,2,3,.…), 


tan( Qi.2h;) = 


ka 2 
tan(g,y hy) = 一 一 一 一 (mz1,2,3,:55;i1 22,3, 
Oy 2 Ga 
Su No -1), 
Qi 26. 2 
Q,;»-arctan| —, — — 
Caio 


其 中 ,B B, s Buw2 为 一 系列 待定 系数 。 
2.2 桩 振动 方程 求解 
将 式 (3) 两 边 同 时 进行 Laplace 变换 ,可 得 
9^ $,(z,s) 
az? 
2mr? oF,s(z,s) =p, Jus Bi(z,s) (26) 
AP, Fii s). Fio (2,8) 和 0; (2,5) 4r HJ 
fia) a (8E) WI in, t) 的 Laplace 变换 式 ,其 中 
fia Gs) Mia) TAERA) X C13) ME. 
1) 当 i=1 时 ,对 桩 的 边界 条 件 式 (12) 式 (13) 
进行 Laplace 变换 ,并 代入 式 (26) ,可 得 


Gd y -2ar; F; (2,5) 一 


2 
2mri ， 


G ` 
G s1,2 
pl Ja 


£ Bm 2d 2l (0,2712) sin[ Q2 (z = h, ) ] (27) 
1 
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式 (27) 的 通 解 为 
Di(z,s) = pie 6y* + Ore are 十 Ys A m aSin[ a4 (z -h,)] + > 7 各 Basin[aw a(z -h)] (28) 
m= m=1 
式 中 将 式 (10)、 式 (11) 进行 Laplace 变换 ,并 将 
yz LTT a Ga agm me (Ina 14) XX (20) 式 (28) 代 入 式 (10), 式 (23) 式 (28) 代 入 
g J ( Caan *tp,s ?) í 式 (11) ,可 得 


2 
(2) _ 2mri 264,38, 2b (Q5 3712) 
ES 2 2 
Ja ( G 05.2 *tp,s ) 


YA. IK quani 1)sinLo,, a(z -h)l= 


m= 


oo oo 
a) NES (2). -, / EE. (1) , (2) : 
T, {a e um * B e M 十 Y Ymi A m SINL a, (z = h, ) ] 十 > Ym B, osinLa, s (z = h, )1 (29a) 
mzl m=1 
EN ali (ami aia) sinLa, s (z - h, )] = = 
m= 


/2 ES Pe 
4 na [Bie ey * Be 三 * Yos ml Sin[ o 1(Z h, )] + EXE ml asino, a B h 而 | 0 


加 假设 cy = Ani = Ani XII (29a) (29b) 积分 ,利用 sin[ qi(z 一 ha)] 的 正 交 性 ,可 得 
两 边 同 时 乘 以 sn[ anı (2 一 hb)], 并 对 z 在 [0,h] 上 


rial (am, 3513) " IT 
Am = EE MES [. (gi"e Ve * BÜ ev Sw) sinLg, (z — h,)1dz (30a) 
PK, /Bw 
B,» = rni : op (B! De a” t Bie T 人 isJsin[a (z xs h, ) ] dz (30b) 


ml ml 
XH 
h1 
La = n sin Lo, (z - h) ldz, E m = ER Cau atii) - y ri MiC aria?) my Ti aK Camara) 
AEA HEE ff A 
Bi(z,s) = Bi? (eV i* + Dénsinto (a-h) 1} BP (eV + Desinta (2 zai 3D 


m=1 
Xm 
r L 5f hl CE K, hl £N 
Ém = YQ A ur E Cpr “sinl a, 16: mx h 1)1dz + yo aa Ear : n, naf Gana) f Cy "sinl a, 16: É h ,)1dz, 


ml Ln E ml 


r4; jr hı rs TiK; r hı Po 
Lm = y9 1,1 u Pe C "sinL o, (z — h, ) dz +y 2) 1,2 iaa aa Gana. E "sinL a, (z — h,) 1dz 


ml 


LE ml Li E, 


结合 式 (8) 和 式 (9) ,得 到 第 1 桩 段 顶部 动 阻抗 为 


(1) ic 
mE ) FUNER 十 È Eni amicos omh n)]* 中 [es 十 DIN j] - 


GJ 


kr = 


(0, "E C plt pl " 
gol - Y Ensin anh )]* [1 一 Y isin lanh; )] 
式 中 BO /B® E MCN 法 ,利用 位 移 连 续 条 件 , 并 结合 sin( Aniz t On) 的 正 
ee NES 可 获得 第 ; REEL HITLER 


qp. (z, s) = ES p” [eV + È EiSin( Qiz 十 Pmi) ] * B® [ : er 十 > Cisin( Qiz 十 ou)] (33) 


AP 
Emi = Ys 
r; L mi, r; hi - / £P . lj K, mi,l' i, ha JE 
É mi = x : (d s ^[ e “i sin( Q miZ + Pmi ) dz + 2 2 Can 4 Ga Mean f'- Cpi wsin( Q miZ F Pmi) dz, 
0 


Yni 


r,l T; 4) ti fev. ER; hi 
D Mua afi uaria) ( e TN z+ pi) dz + yt + er ei et pe debi d: 


Lpi Eu mi ^mi 


E Ls LE Li E, 
2 2 
(dy n. 2mr; 16, is AK (uaria? 2) o 2mr; ne si 2d mi ab(quaria) m M g 
Ymi = 2 ^Ym 二 一 2 QLQ7 sin (Qiz + ni) dz, 
Ja Gua mi,l + p,$ ) Jy Cs + pops ) 0 
Eni = ER Cati) - yC ri Mi(quiari2) -yk T aK, (miari) 


桩 项 动力 阻抗 函数 可 以 表示 为 


— | /Rr +|_ /Pr 
ni Bol p 十 X inset] | e 4 Xue] m 


ks gr EL Ed GJ (1) m 


'" $400, : l 
gal! t X since )] t [1 + » É ,;SIn C, )] 


mii 


EM | p /Bp / Ghi 
C - Cy Jl GC. sev ^i 2 * > Enia Q miCOS( CQ hi; T Omi?) * kal e Z T > £i sint Amih i + Qi )] 
VG m= 


BS 
[ g? m Cu L| =se As F È nanco Ca h; +o) |* kule -yip + Y busina, h; + eu) jy 


m- 


, Ner E i 桩 段 桩 项 阻抗 3) 同 理 , 当 i=N 时 ,利用 cos(Qwz) 的 正 交 性 
OR i+1 桂 段 桂 底 反 力 系数 大 小 相等 , 即 kuoa = 可 得 第 NN 桩 段 桩 身 扭转 角 为 
6 夏 同 时 可 得 第 i+1 桩 段 桩 底 反 力 系数 。 


Bn(z,s) = pi [e ew 十 b» £, cos Co, uz] 2 [e V à V 6 + X daeosta, "m (35) 
m=1 m= 


Ev. 


: hN hN 
AN = t 2ry A (Annana) (^ / rss i o) 27rw2 天 《quw rw a). (^ /名 = i 
LA É mn = Ymy h E e PN cos( Oz) 十 Ymy hE e piy cos( Oz) , 
N mN 0 N” mN 0 


= 27v ， L r hN 2r K, r an 
Ciu = yR ied E xD. EL "cosa, z) dz + y9 — "ow ~, aad EE “cos( a, yz) dz, 


2 
2mry m 1 nv E E arva? (2) 2a sU alos ( Mr 


(1) N,1 YsN,1 dmN, Nin, i 
= , 


Y nN m 2 p, ,»/mN 2 2 
JC G NON ,1 十 PpS ) JC G NON 2 * PpS ) 


三 (1) (2) 
Emn si LK, ( dmnN,IrN,l ) — Ynn Tya M, CA mN oT 2? C YnN ryaK, q.i Ta ) 


桩 顶 动力 阻抗 可 以 表示 为 
(1) p, p, 
my(s) e NGw NG 
一 人 = (36) 


Seda 
BD\0,s) Gon pN 人 
E (1 十 P5 


m- 


k TN 一 


式 中 
o 
= [eaa] - ps Pr cq AN - b £ va, ySin(a, vh j* kanle Vm + X costas Jy 


=a 
Bx 
/ Pp : A TEB 
1| ^ ca a C B Y EunQnvsin( Oc Toy ) | Fw [ e fn * X £,.cosC Am nhy u | 
pN m=1 


m= 
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以 上 为 Laplace 变换 域内 的 解 ,采用 Crump ^ 
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O 
构造 的 Laplace 逆 变 换 公 式 可 得 其 时 域 解 为 


全 RetFcoy 3 [ReF(e + Eti): stw] 三 [mef。 + jsin e] (37) 


式 中 :Ni 为 级 数 截 取 项 数 ;i= V -1;c 与 7 为 计算 参 
ZH Tzu," T-z2ncT =5。 试 算 表 明 , 当 NUR 
50 时 ,能 很 好 满足 精度 要 求 。 


3 数值 计算 与 讨论 


3.1 与 已 有 结果 对 比 


1 与 理论 解 对 比 

管 桩 内 径 趋 于 0， 将 上 述 解 答 退 化 到 实心 桩 
Ms Militano 4&'?' H lits IRD RE ETTXT EG 
m BEIBRERI QE, sS 
JP t -t/G. prii RERNE o = 
m. pr 750, n=0, pzp, 22650 kg/m, k= 
16 m/s, H 25 m, G, ; 2 20 MPa,p, 2 1.3p, n, = 
0X54n,7r, ; 210 ^m, 
co ml 2 可 知 ,本 研究 解 与 Militano 解 描 述 的 规律 
放浪 一 致 ,在 峰值 及 后 期 响应 上 存在 一 定 差异 。 需 
D UE 的 是 , Militano 5x! 的 研究 基于 Novak 薄 层 法 
,不 考虑 土 体位 移 与 应 力 沿 深度 的 连续 性 ,所 以 在 
范围 内 使 得 桩 顶 动力 刚度 偏 小 而 阻尼 人 往 大 。 由 于 
让 信号 实际 上 是 由 不 同 频率 成 分 的 简 谐 振动 结果 蕾 
加 而 成 ， 使 得 薄 层 法 得 到 的 频 域 结果 经 Laplace 逆 变 
换 提 的 时 域 结 果 失真 ,而 本 研究 采用 连续 介质 模型 ， 
可 得 土 体 在 各 频率 下 的 真实 响应 ,结果 更 为 准确 。 


— 本文 解 GVG.=200 
。 Militanoff£ " 

--- 本 文 解 GJ/G=1 000 
* Militanofj£" 


图 2 ZAER AER EG 


Fig.2 Comparison with the results of triangular loading 
3.1.2. 与 数值 模拟 结果 对 比 
图 3 WEHE - 土 相 互 作 用 有 限 元 模型 。 模 型 范 
围 取 60 倍 桩 径 , 桩 周 土 和 桩 芯 土 均 采用 弹性 本 构 模 


型 , 土 体 四 周 为 无 限 元 吸收 边界 , 管 桩 与 土 体 采 用 三 
维 8 结 点 缩减 积分 单元 (C3D8R) ,并 采用 动态 隐 式 
分 析 步 。 桩 土 接触 面 采用 tie 约束 , 土 体 及 桩 端 为 固 
定 约束 。 参 考点 设 在 桩 顶 圆心 处 ,并 与 桩 项 进行 运 
动 耦合 约束 ,在 参考 点 施加 三 角形 , 半 正 弦 及 和 矩形 冲 
击 荷载 。 参 数 取 值 同 下 述 3. 2 节 算 例 ,通过 对 比 数 
值 模拟 及 本 研究 解析 法 得 到 的 桩 项 转角 位 移 时 程 曲 
线 , 发 现 两 者 结果 吻合 度 较 好 ,验证 了 本 研究 理论 推 
导 正 确 性 。 由 图 4 进一步 看 出 ,在 相同 作用 时 间 及 
丛 载 峰 值 条 件 下 ,不 同 荷载 形式 对 桩 顶 扭 转角 的 影 
响 不 尽 相同 ,作用 和 矩 形 荷载 时 的 转角 普遍 较 大 ,三 角 
形 最 小 , 半 正 荷载 介 于 两 者 之 间 。 


^ ' ' 
| Xm o: : 
NES 边界 i 60r. : 
qe 
A i i 
i ! i 
i i i 
[ i 
$ i 

i 


图 3 EIE- 土 相互 作用 有 限 元 模型 


Fig.3 Finite element model of pipe pile-soil interaction 


pi(10 rad) 


tims 


JA 与 数值 模拟 结果 对 比 


Fig.4 Comparison with the numerical results 


3.2 成 层 土 中 管 桩 瞬 态 扭转 振动 特性 


以 下 通过 数值 算 例 分 析 双 层 饱和 土 中 管 桩 瞬 态 
扭转 振动 特性 。 选 取 三 角形 瞬 态 荷载 ,其 Laplace 变 


换 为 :T(s) =27 (1-2 e73)/(n) (to 为 荷载 持续 


时 间 ,io/2 为 最 大 荷载 中 对 应 时 间 ) 。 参 数 取 值 为 : 
土屋 ,H=10 m,N=2,h,/h, = 1, Go / Gaj si kas 
0. 1G, ;,n 20.4, p, 21000 kg/m? ,p, =2 650 kg/m*, 
TE 107 m/s; f ft, T, = 1 000 kN 
rj, =r 20.5 m,r, , 2r, 20.38 m,G, - 10 GPa, 
G,/G, 21000,p, 22 500 kg/m , 且 桩 芯 土 完全 填充 
于 管 桩 。 

为 验证 采用 分 布 式 线 弹 筑 模拟 成 层 土 层 间 动 力 
相互 作用 的 合理 性 和 精度 ,以 单 层 饱和 土 为 例 ,将 其 
划分 为 同性 质 的 双 层 地 基 ( 桩 土 参数 同上 ) ,采用 本 
研究 简化 模型 进行 计算 ,并 与 严格 解 进行 对 比 。 试 
算 发 现 ,不 同 土 层 层 间 分 布 式 弹 得 系数 ,得 到 的 管 
桩 扭转 时 域 响应 结果 与 单 层 均 质 土 的 严格 解 基本 一 
至 一 说 明 ,对 管 桩 扭转 振动 响应 基本 没有 影响 , 简 
化 模型 具有 较 高 的 精度 ,这 也 与 文献 [11] 中 有 关 ,， 
对 联 紧 向 振动 的 结论 相 一 致 。 限 于 篇 幅 , 此 处 未 列 

关 计算 结果 。 
COE 5 反映 了 荷载 作用 时 间 对 桩 顶 及 桩 身 响应 的 
影响 。 由 图 5(a) 可 见 ,增加 荷载 作用 时 间 , 述 顶 所 
& 铺 的 响应 时 间 随 之 增 大 ,其 峰值 逐渐 增 大 但 增幅 
河源 减 小 。 由 于 印 载 后 振动 引起 的 波 仍 可 在 管 桩 中 
人 
和 痊 载 作用 时 间 越 长 ,振荡 越 微弱 。 此 外 ,还 可 看 
出 最 大 桩 顶 所 转角 产生 时 刻 滞后 于 最 大 荷载 对 应 
IA RARR, 当 作用 半 正 弦 荷 载 时 也 存在 上 
述 规律 。 瞬 态 荷载 达到 最 大 时 桩 身 扭转 角 随 深度 分 
布 规律 如 图 5(b) 所 示 。 由 图 可 看 出 ,冲击 时 间 越 
短 ( 桩 身 扭转 角 沿 深度 方向 迅速 衰减 至 零 。 

图 6 和 图 7 为 管 桩 剪 切 模 量 不 变 时 , 桩 土 模 量 
比 对 桩 顶 及 桩 身 响应 的 影响 。 由 图 6(a) 可 看 出 , 增 
大 桩 土 模 量 比 ( 土 体 越 软 ) , 桩 顶 扭转 角 峰 值 随 之 增 
大 ,种 载 后 反应 持续 的 时 间 越 长。 这 是 因为 土质 较 
软 时 , 秃 载 后 波 在 管 树 中 的 传播 不 易 衰 减 , 故 曲线 波 
动 较为 明显 。 同 时 由 图 6(b) 可 知 , 桩 土 模 量 比 越 小 
(土质 越 硬 ) , 桩 身 扭矩 沿 深度 衰减 越 快 ,在 管 桩 上 
半 段 内 就 基本 衰减 至 零 。 桩 土 模 量 比 C, 6, = 
1 000,200 时 ,扭转 角 随 时 间 和 深度 的 变化 如 图 7 所 
示 。 可 以 看 出 , 当 桩 土 模 量 比较 小 , 即 士 质 较 硬 时 ， 
桩 端 附 近 的 扭转 角 已 接近 于 零 ,而 土质 较 软 时 , 桩 身 
扭转 角 变 化 相对 平稳 , 且 桩 端 附近 仍 存在 扭转 角 。 
这 表明 软 土 相 比 于 硬 土 对 桩 的 约束 作用 减弱 , 从 而 
桩 顶 荷载 可 传递 至 较 深 位 置 。 


“mi 220 ms, 
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7.5 
6.0 " X | \! =30 ms 
FA 1 
/ 4 
^ y J \ \ t=20 ms 
4.5 i AJ 
= "TI / \ 
E j |/ \ 10 
3 HOMI E ! ms 
,30F iS ' \ \ 
= it L / ' \ f 
& QUA 1 Y 
EAE 
VN ! 
i E | ， 让 NE i 
;1 dl] ue? er omen : 
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-Eg 
0 4 7 1012 1720 30 30 a 
tims 
(a) 桩 顶 捏 转 角 
0 pa 
a e 
a 


2 di 
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4 P4 A Nt=10ms=5 ms 


1,20 ms,{=10 ms 


/ f 
/ / E " 
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l 0 l 2 3 4 5 6 
pi(10 rad) 
(b) 桩 身 扭 转角 (最 大 荷载 对 应 时 刻 ) 
图 5 和 荷载 作用 时 间 对 桩 项 及 桩 吴 响 应 的 影响 
Fig.5 Effect of duration time of triangular loading on 


dynamic response of pipe pile 


«wl (10 rad) 


"0 s 10 1s 20 25 30 35 40 
t/ms 


(a) 社 项 扭转 角 


GIG, aw 


4 i 
GJG, =1000— \ kd 
6 \ UY. 
N \ I 
A li 
\ i 
8 n 
1-10 ms li 
| | 
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MI(IO kN * m) 
(b) fit & fb 


图 6 桩 土 模 量 比 对 桩 项 及 桩 身 响应 的 影响 


Fig.6 Effect of pile-soil modulus ratio on time domain response 
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v (10 rad) 
ct Oc tX Roo 0 


Nou be 


q/(10 rad) 


O= Oc = 


10 0 
(b) HELER H200 
图 7 桩 土 模 量 比 对 桩 身 扭转 角 的 影响 
Fig.7 Effect of pile-soil modulus ratio 


.00405v1 


on the twist angle along the pile 


CÓ I s Bc TIRHLBEOR n REBEL f B 
Davey 渗透 系数 局 的 变化 规律 。 由 图 可 知 , 在 较 高 和 
镜 攻 ,范围 内 ,相反 所 续 角 随 渗透 系数 几乎 保持 不 
Le, , 较 大 时 桩 顶 扭转 角 大 于 , 较 小 的 情况 ; 
"A E 107 ~1.0 m/s ZEIT, CIA S IRE, 
扭转 角 随 之 增 大 。 这 是 由 于 此 时 土 体 可 看 作 开 
en ad 
的 能 力也 越 弱 , 故 管 桩 扭转 动力 响应 越 明显 。 此 外 ， 
EGTLBIR n 对 应 的 桩 项 扭转 角 也 较 小 , 受 渗 适 系 
影响 也 较 弱 。 


30 


Ios dka) 
3 4 


一 一 /=-0.5 
zer 1-03 Jos a(k,) 


de Ik.) 
t=10 ms 


PIO rad) 


joa) 
8 桩 项 扭转 角 随 渗透 系数 的 变化 规律 


Fig.8 Variation of the twist angle of pile top 


7j 


with permeability coefficient 
这 表明 对 于 密实 土 体 ,孔隙 水 的 参与 作用 将 减 
弱 , 桩 受 土 的 约束 作用 增强 ,从 而 使 得 管 桩 桩 顶 扭转 
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角 减 小 。 由 于 管 桩 在 沉 桩 过 程 中 桩 周 土 向 旁 侧 挤 
压 , 致 使 上 部 桩 蕊 土 形 成 浮 土 ,hi, 减 小 而 wes 
由 图 可 知 , 当 kai # s LESE, ECISEBERS E18 38 2C kaa 
对 桩 项 扭转 角 几 乎 无 影响 。 

双 层 土 剪 切 模 量 差异 对 桩 顶 扭 转角 的 影响 如 图 9 
所 示 。 由 图 可 知 , 增 大 上 层 土 剪 切 模 量 使 最 大 桩 项 扭 
转角 大 幅 减 小 ,外 载 后 桩 顶 反 应 持续 的 时 间 变 短 。 而 
改变 下 层 土 剪 切 模 量 对 桩 顶 响应 的 影响 很 小 ,说 明 管 
桩 瞬 态 扭转 振动 特性 主要 受 上 部 土 层 性 质 的 影响 。 

10 


G,,-G...-10 MPa 


G., 75 MPa, G,.,-10 MPa 


5 Ga =20 MPa, G...—-10 MPa 
r G, =10 MPa, G-,=5 MPa 
6 / e G,,-10 MPa, G..,-20 MPa 


pi(10 rad) 


25 
0 10 20 30 40 
t/ms 


图 9 各 土 层 不 同 剪 切 模 量 对 桩 项 扭转 角 的 影响 
Fig.9 Effect of different shear modulus of soil layers 


on the twist angle atop the pile 
对 于 三 层 饱 和 土地 基 , 中 间 土 层 前 切 模 量 对 桩 
号 瞬 态 扭转 响应 的 影响 如 图 10 所 示 。 


0 E 
Tr d 


vt 10 rad) 
(a) 桩 身 扭转 角 


Tm 
Ue 
Nu 


ü b" .、 8k 


l 
— —6G, JG, 7025, G, =0, NA 
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图 10 软 硬 夹层 土 对 桩 身 响应 的 影响 
Fig. 10 Effect of soft-hard We soil on 


time domain response of pipe pile 


参数 取信 为 : 土 层 总 厚度 女 =10 m, 其 中 下 层 寺 
h =3 m,'hli]- EJ h =3 m, EJ h =4m, 其 他 参 
数 同上 。 可 以 发 现 , 当 夹层 土 为 硬 士 时 , 桩 身 扭转 角 
整体 减 小 , 而 软 土 时 , 桩 身 扭转 角 整 体 增 大 。 由 
图 10(b) 可 知 ,在 土 层 分 界面 (深度 4m 和 7m) 处 ， 
由 于 土 体 剪 切 模 量 突变 ,改变 了 十 体 对 管 桩 的 约束 
作用 ,致使 柱 身 扭矩 分 布 烛 线 出 现 转折 点 。 相 比 于 
均 质 士 , 硬 夹层 使 得 桩 身上 部 扭矩 增 大 ,而 桩 身 下 部 
扭矩 减 小 , 且 扭矩 在 夹层 范围 内 光滑 连续 过 滤 , 软 夹 
层 情况 则 与 上 述 规律 相反 。 经 试 算 发 现 ,车 增 大 夹 
层 土 的 厚度 及 ,其 剪 切 模 量 的 变化 对 管 桩 扭转 响应 
影响 就 越 显著 。 

由 于 现 淡 混 凝 土管 柱 在 沉 柱 过 程 中 可 能 出 现 缩 
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图 11 缩 径 和 扩 径 对 桩 身 响应 的 影响 

Fig. 11 


Effect of radius reduction and expansion on 
dynamic response of pipe pile 
对 于 桩 身 缩 径 , 由 于 上 、 下 变 截面 处 (深度 为 4 m 
和 5 m) 分 别 存 在 相对 缩 径 和 相对 扩 径 ( 即 桩 径 恢 
复 ) ,致使 桩 身 扭转 角 在 上 半 段 和 下 半 段 分 别 增 大 和 
减 小 ,而 桩 喘 扩 径 下 扭转 角 和 扭矩 普遍 较 小 。 因 此 ， 
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相 比 于 桩 身 缩 径 , 扩 径 对 桩 身 响应 的 影响 相对 较 弱 。 
此 外 ,通过 改变 变 截 面 位置 ( 深 度 2m 和 6 m) 试 算 
发 现 , 当 变 截面 越 靠 近 桩 顶 时 , 缩 ( 扩 ) 径 对 桩 顶 转 
角 以 及 变 截 面 处 转角 和 扭矩 的 影响 也 越 大 , 且 随 着 
变 截 面 长 度 的 增加 ,影响 也 逐渐 增 大 。 


4 结 论 


本 研究 利用 分 布 式 线 弹 簧 简化 土 层 层 间 相互 作 
用 ,研究 了 瞬 态 扭转 荷载 作用 下 层 状 饱和 土 中 管 桩 
的 动力 响应 特性 , 获得 了 桩 土 耦合 振动 的 解析 解 。 
所 得 研究 结论 如 下 。 

1) 作 用 和 矩形 荷载 时 转角 值 普遍 较 大 ,三 角形 荷 
载 最 小 , 半 正 弦 荷 载 则 介 于 两 者 之 间 ; 施 加 三 角形 或 
半 正 弦 荷 载 时 ,最 大 桩 顶 扭转 角 产 生 时 刻 均 落 后 于 
最 大 荷载 对 应 时 刻 , 且 最 大 桩 项 扭转 角 随 荷载 作用 
时 间 的 增加 逐渐 增 大 但 增幅 逐渐 减 小 。 

2) 增 大 上 层 土 剪 切 模 量 使 桩 项 扭转 角 峰 值 大 幅 
减 小 ,而 改变 下 层 土 剪 切 模 量 对 桩 项 啊 应 影响 很 小 。 

3) 在 较 高 和 较 低 渗透 系数 范围 内 , 桩 项 扭转 角 
几乎 保持 不 变 , 而 在 中 间 范 围 内 ,扭转 角 随 渗透 系数 
的 增 大 而 增 大 ; 且 较 大 孔 际 率 对 应 的 桩 顶 扭 转角 也 
较 大 ,渗透 系数 影响 越 显著 。 

4) 与 均 质 土 相 比 ,人 硬 夹层 的 存在 使 桩 映 扭转 角 
整体 减 小 , 且 桩 身上 部 扭矩 增 大 ,下 部 扭矩 减 小 , 而 
软 夹 层 下 情况 则 相反 ; 桩 身 缩 径 比 桩 身 扩 径 对 管 桩 
瞬 态 扭转 动力 响应 的 影响 更 大 。 
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